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Abstract 

It has been researched the kinetics of the thermal equilibrium's estab- 
lishment in an early Universe under the assumption of the recovery of 
interaction scaling of elementary particles in range of superhigh energies. 
The case of the thermal equilibrium's weak initial violation and basic cos- 
mological consequences of the thermal equilibrium's violation have been 
researched. 



1 Introduction 

In the previous paper of one of the authors [1] it was shown, that in the case 
of the scaling behavior of the particles' cross-section of interaction in range of 
superhigh energies: 

Const ,^ , 

CTtot ^ , (1) 

s 

where s is a kinematic invariant of the four-particle reaction (details see in 
[T]), the initial particles' distribution in the expanding Universe is not to be 
equilibrium, but can be random. In this paper we investigate the kinetics of the 
processes with elementary particles in the early Universe under the conditions of 
scaling of interactions with the purpose to clarify the boundaries of randomness 
of the initial particles' distribution. As the cross-section of elementary particles' 
interaction at that we will use an asymptotic cross-section of scattering, UACS, 
incorporated in papers [2], [3^1 : 

27r 27r 

f^o(s) = —7 ; r- = (2) 

s (l + In" ^ j ^Ms) 

where sq ~ 4 - the square of the total energy of two colliding Planck masses, 

A(s) = l + \n Const. (3) 

So 



^ As in the previous paper we will use a system of units G = h = c = 1. 
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2 Kinetic equations for superthermal particles 

2.1 A simplification of the relativistic integral of collisions 

The relativistic kinetic equations for homogenous isotropic distributions fait-,p) 
have the form (see the previous paper [T], details see in [i], [5]): 

dfa a dfa 1 ^ ... 

^7 = ^ „ „ 2^ Jab±^cd[t,P), (4) 

dt a dp ^ml+p^ 



where a{t) is a scale factor of the Friedmann's world: 

ds^ = dt^ - a^{t)[dx^ + p^(x)(sin2 9dip^ + de"^)]; (5) 

p' = -5«/3pV, (a,/3=0); (6) 
Jab^cd{t,p) is an integral of four-particle reactions [B], [7]: 

Jab{t,p) = j dnbdn,dnaS^^\Pa + Pb - Pc - Pd)x 



X [(1 ± fa)il ± fb)\fcfdM,d^abV - (7) 



-{l± h){l± fd)fafb\Mab-,cdV], 

characters ± correspond to bosons (+) and fermions (— ), Mi^f are invariant 
amplitudes of scattering (line means the average by particles' states of polar- 
ization), diTa is a normed differential of volume of the a particle's momentum 
space: 

, , — Padp^dp^dp^ , 

d-Ka = V-ff JTT^ , (8) 

Pa is a factor of degeneration. 

Let us simplify an integral of four-particle interactions ([7]) , using the prop- 
erties of the distributions' isotropy fa{t,p). For the fulfilment of two inner 
integrations we will proceed to the local cm. system, integration in which is 
carried out by the elementary way. After the inverse Lorentz transformation 
and conversion to the spherical system of coordinates in the momentum space 
we obtain ([8 ): 

oo 

25*6 + 1 f qdq 



2Sb + l f 

^^'^^ = -W^P J 



sds 





27r 



dt\M{s,t)\^ J dipx 

{fa(p4)fb{q4)[l ± /c(p4 - A)][l ± ,Uq4 + A)]- 
-/e(p4 - A)/d(<Z4 + A)[l ± /a(p4)][l ± Mq4)]}, 
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where 



A2 



/ \ "'a 
P4 - 94 - [Pi + 94) 



ts 



ts 



Apiqi 1 



1/2 



A2 ^ s2 _ 2s(m2 + m2) + (m^ - m^)^ - 
a function of a triangle, s,t are the kinematic invariants (see [T]), 

It is necessary at that to bear in mind the definition of the total cross-section 
of an interaction [S], (see also the previous paper [T]): 
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CTtot 



167rA2 



dtF{s,t), 



where we denoted as is generally accepted: 



Fit,s)^\M{s,t)\^. 
In the ultrarelativistic limit: 

Pi -5 X 2 

> 00 ^ — 7j 00; A s , 

rrii mf 

aforecited expressions are greatly simplified: 



(9) 

(10) 
(11) 



Jab{p) 



{2Sb + 1 

32(27r)4p 



00 ipq 



dq 



ds 



2ir 



dxF{x,s) / d(px 





{fa{p)Mq)[l ± Mp - A)][l ± fd{q + A)]- 
- Mp ~ A)f,{q + A)[l ± UpW ± Mq)]}, 
where a new variable was incorporated (see [I]): 
t 



and 



A = x{p — q) — cos (^■\/ x{l — x){Apq — s). 



(12) 
(13) 
(14) 
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2.2 An integral of collisions for a distributions' weak de- 
viation from the equilibrium 

Let us investigate at first a weak violation of ttie tfiermodynamical equilibrium 
in the hot model, when the main part of particles, no(t), lies in the thermal 
equilibrium state, and only for the minor part of particles, ni{t) the thermal 
equilibrium is violated (see Fig. [T]): 

ni{t) < no{t). (15) 




Figure 1: The schematic representation of the distribution function's deviation 

from the equiUbrium. 

Henceforth in this paper we will suppose that distribution functions different 
not greatly from the equilibrium ones in the range of small values of energy, 
smaller than certain unitary limit, p — pq (or T — Tq), below which scaling is 
absent, and can be violated at energies, above the limit: 



faip) 



fO 

J a 



1 



exp(- 



-Ha + Ea{p) ^ 

T 



P < Po] 



± 1 



(16) 



I A/„(p); /0(p) « Afaip) « 1, p>Po, 



where /Xo(i) are the chemical potentials, T(t) is a temperature of the equilibrium 
component of plasma. Thus, in range p > po it can be observed the anomaly 
great number of particles in comparison with the equilibrium one, but minor at 
that (. (jlSp ) in comparison with the total number of equilibrium particles. 

Let us investigate the process of relaxation of the distribution fa{p) to the 
equilibrium faip)- The problem in such a setting for the special case of the 
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initial distribution f{t — 0,p) was solved earlier in [2], [5]. Here we will give 
the general solution of this problem. At that, as it will be obvious from the 
further arguments, the cosmological plasma can be formally considered as a 
two-component system- equilibrium with the distribution f^{t,p), and nonequi- 
librium superthermal, with the distribution Sfa{t,p) = where number 

of particles in the nonequilibrium component is small, but energy density, as a 
matter of fact, is random. 

Let us investigate an integral of collisions (|12l) in range 

P > Po » T. (17) 

In consequence of the inequality (jl6p in this range we can neglect collisions of 
superthermal particles between themselves, not going beyond the account of the 
superthermal particles' scattering on equilibrium particles. Therefore the value 
of one of momentums in the integral of collisions, P' = p — A, or q' = q + A 
must lie in a thermal range, another's value - in a superthermal one, outside 
the unitary limit. A subintegral value of the integral of collisions is extremely 
small outside this range. In consequence of this circumstance we can neglect 
the second member in curly brackets (jl2p . since it can compete with the first in 
asymptotically small variation ranges of variables x and ip: x{l — x) < T/p 0. 
The statistical factors of type [1 ± fa{p')] in the first member of an integral (|12p 
can noticeably differ from one besides only in the range of momentums' thermal 
values. As a result the integral of collisions p2)) in the investigated range of 
momentums' values can be written down in a form; 



(25b + l)A/ a(p) , 
{2ttYp 



'Jab-ir^cd[P)\p-^pg — /r,_\.ci„ ^ 



oo 2p{q^+q) 1 

^fim^ f ^ dxFM (18) 
J Jmf+cp J IGtt J 

^ ° 2p{qi-q) 

Using the definition of the total cross-section of scattering (|12l) , we obtain from 

Jabi^cd{p)\p^p^ = 

(2g, + l)A/.(p) f qf^{q)dq f 
= TTT^ / / o / <ytots{s)ds. (19) 

^ 2p{q^-q) 

Substituting eventually an expression for utot in the inner integral in form of 
UACS, ([2]), carrying out an integration with a logarithmic accuracy and sum- 
ming up the obtained expression by all channels of reactions, we find finally: 

'^a(p)lp>po ~ 

oo 

.-A/.(p)Ve^t±iliii/4^^, (20) 
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TT 



^mi + q2 A(s) ' 
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where 

- 1 4 

Vab is a number of channels of reactions, in which a sort particles can participate 
a. 

Let us calculate values of the integral (I^Hl) in extreme cases. 
A scattering on nonrelativistic particles. If h sort equilibrium particles 
are nonrelativistic, i.e., q <^ mt, the integral (j20l) is reduced to the expression: 

=-s.'aup)T.^ , (-^>^)- (21) 

b 1 + In 

A scattering on ultrarelativistic particles. If b sort equilibrium particles 
are ultrarelativistic, i.e., mf, <^ T, moreover their chemical potential is small, 
- Hb ^ T, then calculating the integral (l20l) with respect to the equilibrium 
distribution (|16p. we find: 

= '3 i + in2y^^^ A/4p). (^f « Mb « T), (22) 

where 

^(25+l) + i^(25+l) 
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Nb is a number of sorts of equilibrium bosons, F - fermions. 

To estimate contributions of nonrelativstic and relativistic particles to the 
integral of collisions in an equilibrium component, we first will calculate their 
concentrations. A concentration of ultrarelativistic particles in the hot model is 
resulted from the expression (1161) for the distribution function of an equilibrium 
component by means of the substitution 

E{p)=p; tia = (23) 

into the formula for the determination of particles' number density (see [T!): 

oo 

na{t) = J fa{t,p)p'dp. (24) 



Thus we find (see [T]): 

na{t) = V^^' g"C(3) (25) 

A concentration of non- relativistic equilibrium particles at under the assumption 
of conservation of their number is violated in proportion to a^^{t). Therefore in 
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conditions of equilibrium's weak violation the ratio of nonrelativistic particles' 
density to density of relict photons is approximately constant (since T ^ a{t)~^: 



nj{t) 



Const = (5 - 10"'" ^ lO"". (26) 



Calculating the ratio of contributions to the integral of nonrelativistic and ul- 
trarelativistic particles' collisions, we obtain: 

Jnon Jultra ^ = C{3)S ^ ^ IQ-^ — , 27 

- ratio of contributions is small at T ^ lO^mf, and diminishes with time. There- 
fore in future we will neglect the contribution of nonrelativistic particles in the 
integral of collisions. 



3 A Relaxation Of A Superthermal Component 
On Equilibrium Particles 

Substituting the received expression (|22]) for the integral of collisions into the 
kinetic equations ([4]), we obtain a kinetic equation, which describes an evolution 
of an ultrarelativistic superthermal component in the equilibrium cosmological 
plasma: 



Taking into account the fact, that the variable: 

V - a{t)p, (29) 

is an integral of motion [1], we proceed to variables t,!^ in the equation (pS)) : 
for any function ^(i,p) at that the following relation takes place: 

At such a substitution the kinetic equation (|28p can be easily integrated in 
quadratures. For convenience in future we will define more exactly the normal- 
ization of the variable V . At that there appears a necessity to compare values, 
which are used in the nonequlibrium model, with corresponding values of the 
standard cosmological scenario, since all observed cosmological parameters are 
interpreted in SCS terms. We also should keep in mind two synchronous models 
of Universe: the real - nonequilibrium model M with macroscopic parameters 
P{t) and the ideal - equilibrium model Ai^, which in given point of time t 
possesses certain macroscopic parameters Po{t). 
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An Ultrarelativistic Universe 

Let us consider the Universe with the ultrarelativistic equation of stat(3 (a 
characteristic of a barotropic line is p — 1/3): 

e = Zp. (31) 

Then according to Einstein's equations an energy density of the Universe is 
changed by law: 

ea^ = Const; e = -^^^^ (32) 

and a scale factor is changed by law: 

a(t) - t^'^. (33) 

From the other hand, an energy density of an equilibrium plasma is determined 
via its temperature by the relation (see [T]): 

£o=AA^. (34) 

Therefore, if the Universe was filled up only by equilibrium plasma, its temper- 
ature ro(i) would change by law (see [T]): 

A/-^/^?^o(i)=(^)'^^-^/^ (35) 

- here we take into account a possible weak dependence of an effective number 
of equilibrium types of particles from time, M{t). So, let us define more exactly 
the formula (|29|) by following way: 

P^VN^'^To{t). (36) 

Such is the meaning of the momentum variable V according to this formula: 
to within the numerical factor of order of one V is the relation of particles ' 
energy to their average energy in the same point of time in the locally equilibrium 
ultrarelativistic Universe. 

Thus, solving the kinetic equation (1^51) with the account of relations pop 
and (l36l). we find its solution: 



A/,(t,P) = A/0(P)exp 
where: 



m.V) }yHt')dt' 



V J ^ 



A/0(7') = A/,(0,7'), 



(37) 



It should be noted, that pressure and momentum have the same denotations. 
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is an initial deviation from the equilibrium and there is incorporated the dimen- 
sionless function: 

y{t) = ^ (38) 

and the parameter, weakly dependent from variables t^V: 

A{x) 1 + In^ X. (40) 

Values V ^ 1 correspond to the approximation px ^ T{t). 

Since T{t) is a temperature of the equilibrium component of plasma, and 
To(<:) is a temperature of the completely equilibrium in a given point of time 
Universe, the following condition is always fulfilled: 

y{t) ^ 1. (41) 

To be the correct solution of the kinetic equations, the function Afa{t,V) 
has to satisfy in all times to the integral condition P3)) . Since according to 
the solution (pS)) the distribution function's Afa{t,P) deviation from the equi- 
librium strictly diminishes with time, for the validity of the solution (1361) it is 
sufficient to the function Afa{t,V) to satisfy the condition ((T5|) in the initial 
point of time. This gives: 



oo 

/ 



Af°iV)V^dV » j^yl (42) 



where yo = y(0) < 1. 

As an example we will consider a relaxation of a superthermal component at 
the initial distribution in form of a staircase function for the density of particles' 
number: 

n^AN 

Af{V) = { ^oV^ , (43) 

0, r Vo; 



so that: 



AN^^ I Af{V)V^dV (44) 



is an initial conformal density of the non-equilibrium particles' number. On Fig. 
[2] an evolution of the superthermal component for such a distribution, at that 
we laid y(t) = 1 is shown . 
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Figure 2: A relaxation of a superthermal component for the distribution i[43l) 
with the assumption y{t) = latVo = 100, M /\/N — 10. A relative magnitude 
of the distribution function of particles' number density by conformal energies 
V is shown. Top-down (along the figure's left border) there are firm lines: 
t = 0,t^O,01,t^O,l,t = l,t = lO andt = 100000; a dotted line is t ^ 3. 
Time is measured in seconds. 

Let us remind that cosmological time t is measured in Planck units. There- 
fore a native question, if methods of classical (nonquantized) kinetics can be 
used in times of order of several Planck times, appears. A condition of apphca- 
bihty of a particles' semi-classical description in the cosmological situation is a 
relation, resulting from the Heisenberg's uncertainty relation: 

M > 1. (45) 

According to (|35l) and (|36|) : 

Therefore a condition of applicability of a particles' semi-classical description 
takes form: 

tV'»^|^»2,Q^. (47) 

Such consideration at t ^ 1 is justified for sufficiently great values of the con- 
formal momentum V ^ 1, which exactly correspond superthermal particles. 
Thus, a semi-classical description of particles is applicable in Planck times of 
an evolution of the Universe with the more validity the lesser is a ratio of a 
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thermal energy to a particle's energy. Thus, a description of a superthermal 
(non-equilibrium) component's evolution does not require a quantum consider- 
ation. 



Non-equilibrium Universe 

Let us consider now the Universe on non-relativistic stage of expansion. In this 
case an equation of state is 



p = 

and a scale factor is changed by law: 
a{t) ~ t^/\ 



(48) 



(49) 



We will consider at that a scattering of superthermal particles on equilibrium 
massless relic particles, a number of which in SCS is greater than a number 
of non-relic particles approximately in 10®. Let to is a moment of an ultrarel- 
ativistic equation of state's replacement by a non-relativistic one in the non- 
equilibrium Universe and T^{t) is a temperature of relic photons in the equilib- 
rium Universe. Then: 



45 



1/4 1/6 

[a 

t2/3- 



Integrating the kinetic equation 



we find in this case: 
t 



x{t,V) f y\t')dt' 



to 



A/,(i,7')=A/,(to,7')exp 
where Afa{to,V) is determined by the solution ([57)) . and 

xit,v) 



(50) 



(51) 



7rAAr°(to)tg/' 



AC^VTOit)) 

is a slowly changing parameter, Af is determined for massless relic particles 



(52) 



4 A Heating Of An Equilibrium Component By 
Superthermal Particles 

4.1 An Energy-balance Equation 

In spite of a small in comparison with an equilibrium number of non-equilibrium 
particles, energy, included in the non-equilibrium tail, can turn to be sufficiently 
greater than an energy of an equilibrium component, if sufficient distortions of 
the distribution lie in a superthermal range, to which great values of the mo- 
mentum variable V ^ 1 correspond. Superthermal particles, colliding with 
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equilibrium ones, transfer theirs energy to equilibrium particles and thereby 
heat up an equilibrium component of plasma. For the finding of the real tem- 
perature, T{t), of an ultrarelativistic plasma with the account of its heating by 
superthermal particles we will use the equation (j32[) . which presents itself a con- 
sequence of the energy's conservation law and determines the dependence of the 
ultrarelativistic Universe's energy density, s, from the cosmological time. This 
energy density is composed from the energy density of an equilibrium plasma, 
eo P3p . and the energy density of a superthermal component, Si: 



where 



So + El, 



£1 = + 1)7 ^fa{t,v)r'dv. 



(53) 



(54) 



For the further it is convenient to incorporate a dimensionless variable cr{t): 

£0 



a(t) = ^ = — ^ ^ 1. 

£ £0 + £l 



(55) 



Since the total energy density e{t) is determined from the other hand via fourth 
degree of temperature, To(i), of the completely equilibrium Universe, and an 
energy density of an equilibrium component £0 (t) is determined via fourth degree 
of temperature T{t), an incorporated dimensionless variable is concerned by a 
simple relation with a dimensionless variable y(t) = T(i)/To(i), incorporated 
earlier ([38]): 

a{t)^y\t), (56) 

from which right away the inequality (|4T|) follows. 

Thus from ([53]) subject to the solution ([37| for the non-equilibrium distri- 
bution function we obtain an integral equation relative to function y(t): 



00 



X exp 



at,r) fv\t')dt' 



V J Vt^ 





dv ^ 1. 



(57) 



From this equation in zero point of time we obtain the relation: 

00 

^ J2{2Sa + l)J r'Af,iV)dV = 1 - ao. 



(58) 
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4.2 A Solution Of An Energy-balance Equation 

At specified functions Af^{V) an equation (1571) is always integrated with a 
logarithmic accuracy in quadratures. Actually, instead of variable t and function 
y(t) we incorporate new dimensionless variable r: 

r = XVt (59) 
Vo 

and function Z{t): 

T 

Z{t) = 2 [ v\T')dT', (60) 





where 1 = ^(T'o), 

CXD 

E(25a + 1) / dVV'Afm 

-Po = ^_ (61) 

E(25a + 1) / dVV^AfO{V) 

a 

is an average value of the momentum variable V in the point of time t = 0. At 
that we obtain: 



y=\llz^; Z(0) = 0; Z;(0) = 2V^, (62) 

where uo = 17(0) — y^{Q)- Then after an integral's calculation with a logarithmic 
accuracy the equation ([57)) subject to the relation ((55|) is reduced to the form 



= 2^1-{l-aoMZ), (63) 
where function ^{Z) is incorporated: 

HZ) = m = 

00 

E(2^a + 1) / dPP3^/r(7')e-^^°/^ 
= • (64) 

E{2Sa + 1) / dvr^Af2{v) 

a 

Henceforth it will be convenient to proceed to a new dimensionless momentum 
variable: 
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such that: 



^ = - h ^Po^l- (66) 

E(25a + l)/dpp2A/o(p) 



Then: 



E(25a + 1) / dpp'A/0(p)e 



-Z/p 







E(25a + l)/dppW°(p) 



(67) 



It is obvious, that ^{Z) is a monotone decreasing function of Z, since it is 
always fair 

$'2<0, Ze(0,oo), (68) 
at that <i>(0) = 1 and $(oo) = 0. From this it follows: 

s$ $(Z) ^ 1. (69) 
From the definition (|67|) it also follows, that the following relation is always fair 

<i>^2>0, (Ze(0,c^)), (70) 

therefore graph of function $(Z) is concave. In consequence of a strict mono- 
tonicity and continuous differentiability of function $(Z) according to the equa- 
tion (j66|) function Z[t) with its first derivative are monotone increasing func- 
tions of the variable x. Integrating the equation (|66p . we obtain: 



'\ (71) 

VI - (1 - fTo)<f ((7) 

Stated monotonicity properties of functions ^{Z) and Z{t) guarantee us the 
equation (j7ip at each specified value r has the unique solution Z{t). At a 
computed value of function Z{t) a value of temperature T{t) (or y(t)) of an 
ec^uilibrium component can be obtained according to (|65p and (|66p by means of 
the formula: 

y = [1 - (1 - aomZ)Y". (72) 

Thus, the problem of heating of an equilibrium component of plasma is formally 
solved. 

It should be noted, that a condition of approximation's applicability of the 
LTE's weak violation (|42l) in terms of incorporated here values Vq and (Tq takes 
the form: 

N^'^al'^Va » 1 - ^0. (73) 
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4.3 An analysis of a thermal heating's process of estab- 
lishment 

Let us proceed now to the analysis of the obtained solution. 
An asymptotic behavior at small times 

First we will consider an asymptotic behavior of solutions at small cosmological 
times: 

t ^ =^ T < 1. (74) 



Expanding a subintegral expression in the right side of (|67l) in Taylor's series 
by the degrees of smallness of we obtain subject to the definition ([M)) an 
asymptotic expansion of function $(^): 

$(Z) = 1-Z + 02(Z). (75) 

Substituting ([75]) into the equation (1551) and integrating the obtained equation 
subject to initial conditions (j62p . we find an asymptotic, at small times, solution: 

^-tV^ + t2(1-cto), (76) 
from which subject to ((62|) we obtain: 



The case ctq ^ 1 corresponds to a small energy density of an equilibrium com- 
ponent in comparison with an energy density of a non-equilibrium component 
of plasma in the point of time t = 0. According to ([75)) a consideration of this 
case is justified for sufficiently great values of T'q: 

^0 » (78) 

In this case we obtain from (|77p an asymptotic law of temperature's variation 
in the early Universe: 

T(t)^T^^(t){^^'\'l'^r^l^ (79) 

a temperature of plasma falls more slowly than in SCS, but in each point of 
time the real temperature is lower than the corresponding temperature in SCS: 

y{t) ^ 1 ^ T{t) ^ To(t), (80) 

and is compared with the last in a certain point of time 1. Let us find this 
time. Substituting (|79p into the solution ([57]) . we find a law of evolution of the 
superthermal particles' distribution in the Universe's early stage: 

Ma{t.V) = A/°(7')exp (-|^)- (81) 
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According to ([7^ and (|8ip LTE as a whole is recovered in times 



^>i=(y) , (82) 

is exactly a time of an equilibrium's establishment. A plasma's temperature 
at that reaches its equilibrium value To(i) (constantly decreasing at that with 
time). However, for particles with energies higher than the average one Vo, the 
equilibrium is not reached yet in this time. 

An asymptotic behavior at great times in the ultrarelativistic Uni- 
verse 

Supposing t >t and, consequently, y{t) « 1, we will obtain an asymptotic form 
((37|) on great times: 



(83) 



thus, thermal heating's establishment time for particles with the momentum V 
is: 



t-p 



n 2 / ^ X 2 



Let us suppose, that in range of great values of V the initial distribution of 
superthermal particles is extrapolated via power law: 

Af°iV)-^V-^; (7'»1;A>4). (85) 

Then an average energy of superthermal particles Ei(t), for which a thermal 
equilibrium to moment of time t is not reached yet, is equal: 

E,it) = at,V)To{t)Vt « ^. . ^ Const (86) 

practically does not depend on time. In the previous paper [1] we noted that 
in the modern stage 1/A ~ a^, where a is a fine structure constant. Thus, we 
obtain from (j86p the estimation: 

Ei{t)r^ IQ-^ (87) 

is in ordinary units Ei{t) ^ 10^^ Gev - i.e, particles with an energy of order 
of a grand unification character energy or higher stay non-equilibrium. This 
conclusion is fundamentally important for cosmological scenarios. 
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An asymptotic behavior at great times in the non-relativistic Universe 



Let us consider now an evolution of superthermal ultrarelativistic particles in 
the non-relativistic Universe, the thermal equilibrium in which on average is 
already recovered. Supposing in this case y = 1 in the expression ([5T|) for the 
distribution function of superthermal particles, we obtain: 



Afait,r) = Afa{to,V)x 



X exp ■ 



V 



1/3 



(88) 



the distribution of superthermal particles evolves more slowly, than (|83|) . A time 
of the establishment of the thermal equilibrium for particles with momentum V 
in this case is equal: 



1 



and an average energy os superthermal particles falls slowly with time: 



Eiit) 



10 



(89) 



(90) 



1/3 



Gev. 



In the modern age this value is of order 10^^ lO^'^Gev at variation of to in 
limits 10^ IQi^gec^ 



Numerical model 

Let us proceed now from estimations to numerical calculations. Let us present 
an initial deviation of the distribution function from the equilibrium in the form: 

^/"(^) = ^3(fc4^x(l--), fc^O, (91) 

where xi^) is a Heaviside's function (a staircase function), x — V/Vq is a dimen- 
sionless momentum variable, A, Vq and k are certain parameters; a parameter 
k is introduced for the security of convergence of all distribution function's mo- 
ments in the range of momentum's small values. The distribution function of 
superthermal particles' energy density at that has a form, similar to the spec- 
trum of so-called flat noise, when all values of energy are equiprobable (see Fig. 
[3] and [10]). 
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Figure 3: An initial deviation of the distribution function of the superthermal 
particles' energy A/°(x)x'^ from the equilibrium by the formula i91]) . 
bottom-up is: A: = 0, 1, fc = 0, 01, fc = 0, 001. 

Calculating a number density of particles relatively to the distribution ([9T|l , 
we find according to : 

711 (i) - 

_ A(2^. + 1)AA3/^ / k 1 

Thus, we obtain in the extreme case: 

To Win- (fc^O). 



ni{t) ~ A- 



from which follows, that for a fulfilment of the condition (|15p in all times it is 
necessary and sufficiently the following relation to be fulfilled: 

A^^/Hnl^2a'J\{3). 
Calculating now an energy density relatively to the distribution fj9ip . we obtain: 

£i{t) = A 5 7^0(1 - k) . 

Thus, subject to p6p we obtain a dependence of a momentum variable's average 
value through parameters Vq and k (Fig. [l]): 



{l-kf 



In- 



VTTP - 1 \/TTP 
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Figure 4: A dependence of variable's relative average Vq/Vq from the 
parameter k of the distribution Ii91\) . 

Then an energy-balance equation ([55)1 or ([57)1 at < — > leads to the relation 
between the distribution parameters (j9ip : 



^A{2S+l)Vo{l-kf = l-aa. (92) 

TT 

Thus, the following condition, imposing a restriction on model's parameters, 
has to be fulfilled: 

^T'o < ^, (93) 

Calculating the function ^{Z) relatively to the distribution ((9T|) according to 
the formula (I57|) . we find, proceeding to the limit A: — > 0: 

+Ei(-Z), (94) 

where Ei(z) is an integral exponential function [11) : 

z 

r e* 

Ei(^) = / jdt, |arg(-z)| < tt. (95) 



The equation ((7T|) accounting the relation ([55)) takes the form: 



^ _ 1 /■ du 



2 7 - (1 - (To)(cxp(-u) + uEi(-u)) ' 



(96) 
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where r — x/Vq. If we could integrate ([Ml) and found the dependence 9{t) — 
Z{x), we could receive the required dependence y{x) from the relations (l62l) - ((63 
in the form: 



y 



l-(l-ao)(e-^(^)+f?(r)Ei(-%)) 



1/4 



(97) 



However, the dependence 0{t) of course, can not be found. Therefore we act in 
the following way: we find the dependence t{9) from the equation (j96p by means 
of the direct numerical integration. Then the equation (j97p can be considered 
as an equation of type y — y{9), and the set of two equations (l96l) and ([QT]) - as 
parametric equations of the graph y = y{T): 



r{0) 



du 



27 v^l - (1 - CTo)(exp(-w) + uEi(-u)) 



(98) 



y{e) = [1 - (1 - a^){e-'^'^ + 0(r)Ei(-0(r))' 



1/4 



On Fig. |6]the results of a numerical modelling of the equilibrium component's 
heating process by superthermal particles for the initial distribution (PT|) on the 
basis of equations ([M)) are shown. 

From this figure it is obvious that at t 1 temperature of plasma's equi- 
librium component in fact becomes saturated up to the value ro(i)j i-C-i to the 
2 

point of time t ^ Vq /£, , that coincides with the asymptotic estimation 
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Figure 5: A relaxation of plasma's temperature to the equilibrium: 
y ~ T{t)/To{t) subject to the parameter ao: bottom-up 
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(To = 0, 01; 0, 1; 0, 3; 0, 5; 0, 7. Vaiues of the dimensionless time variable r are 
positioned on the abscissa axis. 



5 Basic regularities of a cosmological model with 
an initially weakly violated thermal equilib- 
rium 

5.1 Concentrations of ultrarelativistic relic particles 

Equilibrium concentrations of ultrarelativistic particles in the Universe with a 
weakly violated thermodynamical equilibrium are determined by the relation 
which subject to the definition (j3Hl) can be rewritten in terms of the in- 
corporated function y(t): 

na{t) = ^^ff„C(3) = nl{t)yHt) ^ nl{t), (99) 

TT 

where nJJ(t) are the equilibrium concentrations of the same particles in SCS: 

nl{t) = ^g„C(3) = 



p f 45 



2 ' 



3/4 



= ^.9nm[^] t-^". (100) 



Thus, the number of ultrarelativistic particles lying in the equilibrium in each 
given point of time t is lesser than in SCS. If the interaction between particles 
is described by the scaling cross-section of type ([1]) or ([2]), than a number of 
non-equilibrium particles is described by the distribution function of type (P7)) . 
Let us consider particles' reactions in the range of low energies, in which scaling 
can be violated. Let T^j^{t) is an effective time of interactions of sort "a" 
particles with the other particles and let in the investigated range of energies 
the inequality: 

<fi{t)>t (101) 

has as its solution: 

t > tl. (102) 

Then in the point of time t = sort "a" particles cease to interact with the 
other, i.e., become relic — the number of such stable particles in times, later 
than is conserved, but their density evolves henceforth by law: 



na{t)^nl{t:)y'\t:) 



ait) 



(103) 



Therefore, if in the point of time t = t* an equilibrium is not recovered as 
a whole, the number of relic "a" sort particles is in y^{t*) smaller, than it 
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is received in the standard scenaricl^. Thus, varying a parameter of the non- 
equihbrium Universe's model Vo at co ^ we can regulate the number of relic 
particles and make this number an arbitrarily small [5] - at a increase of Vg the 
number of relic particles decreases. Thus, at Pq > 10^ and ctq ^ 1 relic extra- 
massive bosons disappear, at T'o > 3 • 10^ and ao ^ 1 relic neutrinos disappear. 
This approximate estimator, made in papers, [3], [5], will be specified below. 



5.2 Relic Neutrinos 

So, let us investigate the problem of relic neutrinos' outcome during the harden- 
ing process. A thermal equilibrium of electron and muon neutrino is established, 
generally, by reactions: 



e + + l^e, 

A*^ ^ e+ + 17^ + I^e 



(104) 



The cross-section of an clcctro-weak interaction, corresponding to reactions of 
the neutrino annihilation (|104p in the interesting for us range of sufficiently low 
energies i?, is described by the expression: 

^""""W' ^^^^^ 

where w GnU := 1,4358 • 10^'*^erg/cm^ is a constant of an electro-weak 
interaction. Calculating the hardening time of electron neutrinos, t^, by means 
of the formula: 

Teff = / ^ = t„, (106) 

where it is necessary to substitute E — T{t), an expression for an equilibrium 
density of ultrarelativistic electrons (j25|) and a temperature of an equilibrium 
component T{t), calculated above. Then the ratio of the electron neutrinos' 
number after the hardening in a non-equilibrium model to the same number in 
an equilibrium one N^, is determined by means of the expression: 

On Fig. [6] the results of numerical calculations of electron neutrinos' outcome 
in a weakly-equilibrium model of the Universeare are shown. These results, in 
general, confirm referred above qualitative assessments of previous papers. 



^The time t* is often called the hardening time (see for example [12)'1. 
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Figure 6: An outcome of electron relic neutrinos after the hardening relative to 
a standard model, n^/n1, (an ordinate axis) subject to a parameter of a 
non-equilibrium model, IgVo, (an abscissa axis). Top-down: ctq = 0, 9, 

(To = 0,5, (To 0, 1. 



5.3 A hardening of neutrinos and cosmological helium's 
generation 

A thermal equilibrium of non-relativistic neutrinos is established in reactions of 
type: 

=+ _ 

(108) 



1^ -\- n 



p + v, 
p + 



A concentration ratio of relic neutrinos, Nn, is determined via the kinetic equa- 
tioifl: 



dt 



+ (a + h)Nn = b, 



where: 



Nn{t) + Np{t) = 1, 



(109) 



(110) 



(Np is a concentration ratio of protons), coefficients a{t) and b{t) (velocities of 
reactions (llOSp ) are related via the identity: 



6(t) =a(t)e-^™='/^(*), 



(111) 



Am — m.n — nir, 



1, 3Mev is a mass defect of neutron. 



^Details see, for example, in |12| , 
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Let us investigate the equation (|109p . At f — > according to piip a{t) 
b{t), therefore at i -> from (|105p we obtain the relation: 



dN„ 



dt 



a(<)[l-27V„(0)]. (112) 



From (|112|) right away follows, that at Nn{0) > 1/2 the concentration of neu- 
trons from the very beginning decreases with time, and at Nn{0) < 1/2 - in- 
creases with time. In the second case the function Nn{t) always has a maximum 
in the point of time i*, determinable by the relation: 

- ^^^AL/nu) < I- (113) 

At Nn{0) = 1/2 developing into the Taylor's series Nn{t) by the small t and 
using at that the equation (|109p and limitary relations at t — >■ 0, we obtain: 



, , 1 1 , , Atoc^ dT 
Nnit)^- + - ait) 



2 A ^ ' dt 



t^ 



Since the second member in the right side of the development is the negative 
one, then at iV„(0) = 1/2 the concentration of neutrons decreases. 

Let us proceed now to the numerical modelling. The coefficient a(t) in units 
sec~^ at high temperatures is described by an approximate expression (see for 
example, ['12 J: 

a(t) w 1,61 • T'^(t), (114) 

where T{t) is a temperature in Mev. At low temperatures a{t) — W 10^^ sec~^ , 
where is a probability of free neutron's decay. Supposing in (j35p subsequent 
to [12] = 4, 5, we obtain for To{t) the following relation: 

To = 0,89^^1/^ (115) 

where time is measured in seconds. How the numerical calculations show, this 
point of time is reached sufficiently fastH, - the further history of the nucleosyn- 
thesis does not depend in fact from the initial concentration ratio of neutrons. 

Let us consider reactions (|108p on conditions, when the equilibrium in plasma 
is not reached as a whole, i.e., y{t) < 1. Supposing at that function y{t) obey the 
asymptotic form (|77l) . and at interesting for us time scales A ~ ~ (137)^, 
we obtain: 



f «0,89^y 7^ + 0,523(1 -ao)^, (116) 
where: 

P ^" 

^0 



1017- 

On Fig. [TlfTUlthe results of the numerical integration of the equation pi6p are 
shown. 



^The value of Nn at that in maximum is amount to approximately 0,45. 
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Figure 7: An influence of tlie neutrons' initial concentration on nucleosynthesis 
in the none- equilibrium Universe at parameters' value Vq — 10^^ and ctq = 0, 1: 
along the abscissa axis time in seconds is placed, along the ordinate axis - the 
concentration of neutrons, Nn{t) is placed. The thin line is iV„(0) = 0, heavy 
line is Nn{0) = 0, 5, dotted line is Nn{0) = 1; the value Nn = 0, 165 is noted 
via the dotted straight line. How it is seen from the figure, all three curves 
practically coincide, differing only for small values of time, t < 0, 2sec. 

The results of the numerical integration, shown on Fig. [71 confirni the fact, 
that the final concentration of neutrons practically does not depend from their 
initial concentrations. The number of neutrons after the hardening at Vq ^ 10^^ 
practically coincide with the similar number in the hot model (SCS). The final 
outcome of neutrons after the hardening at Vq > 10^7 increases with the growth 
of Vq, however at the constant value of Vq slowly decreases with the growth of 
the equilibrium parameter (Tq- 

Adduced results of numerical calculations of the neutrons' hardening in the 
non-equilibrium Universe, confirm, in general, estimations, obtained earlier in 
papers [3] and [8], merely specifying their in details. Let us note, that a maximal 
weight concentration of the cosmological He4, Nna, is equal to (see [12j): 

max(Arffe4) =2A^n- (117) 
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Figure 8: An influence of tlie parameter Vq on tie nucleosynthesis in the 
non-equilibrium Universe subject to the nucleosynthesis in the non- equilibrium 
Universe at Nn{0) =0,5 and ctq = 0, 1: on the abscissa axis time in seconds is 
placed, on the ordinate axis the concentration of neutrons, Nn{t) is placed. 
The lower line is a classic result (see fT2Fj, the middle line is Vq = 10^^, the 
heavy line is Vq = 10^*; by means of the dotted line the value Nn = 0, 165 is 

noted. 
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Figure 9: An influence of the parameter Vq on the nucleosynthesis in the 
non-equilibrium Universe at -/V„(0) = 0, 5 and do = 0, 9: on the abscissa axis 
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time in seconds is placed, on the ordinate axis the concentration of neutrons, 
Nn{t) is placed. The lower line is a classic result (see Ji2f ). the middle line is 
Vo — 10^'', the heavy line is Vo = 10^^; by means of the dotted line the value 

Nn = 0, 165 is noted. 
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Figure 10: An influence of the parameter ao on the nucleosynthesis in the 
non-equilibrium Universe at iV„(0) = 0,5 and Vo = 10^*.- on the abscissa axis 
time in seconds is placed, on the ordinate axis the concentration of neutrons, 
Nn{t) is placed. Top-down: ctq = 0, 1, ao — 0, 5, ao — 0, 9, ao — 0, 99; by 
means of the dotted line the value Nn = 0, 165 is noted. Last two lines are 
practically indistinguishable. 

5.4 An observation restrictions on parameters of the non- 
equilibrium distribution 

Owing to the concentration growth of the cosmological Ha an observation 
data about the content of He4 in the Universe on the level not higher than 
35% at increase of the parameter Vo imposes restrictions on upper values of 
the non-equihbrium model's parameter Vo at a specified value of its another 
parameter ctq . This fact was noted earlier in papers |3] , [5] , where on the basis of 
quantitative assessments the range of allowed values of non-equilibrium model's 
parameters was specified. However, these values, correct at not large values of 
the parameter ao, to be improved in the range ao ^ 1. Such improved picture 
is represented on Fig. [11] 
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Figure 11: An allowed by the concentration of He^ range of 
weakly-equilibrium model's parameters IgVo (an ordinate axis) ctq (an 
abscissa axis) - an excluded range of parameters' values is hatched. 

Conclusion 

Summing up the paper's results, we note that modern data of an observation 
cosmology do not contradict the supposition about a weak initial violation of 
the Universe's thermal equilibrium in terms of the condition (I15|) . i.e., to the 
initial smallness of superthermal particles' number in comparison with the num- 
ber of thermal ones (see the condition ()15|) . An evolution of the non-equilibrium 
Universe in the case of the initial distribution's weak deviation from the thermal 
equilibrium, is managed generally by two dimensionless parameters: a relative 
contribution to the equilibrium component's energy density ^ (Tq ^ 1 and a 
ratio of the superthermal particles' average energy to the temperature of the 
equilibrium Universe 1 < Vo < +oo in the same point of time. One could say 
also in such a way: the parameter Vo is equal to the non-equilibrium particles' 
average energy in Planck scales of energy on Planck point of time. An obser- 
vation data at that, generally the prevalence of the cosmological -ffe4, imposes 
certain restrictions on the admitted range of weakly-equilibriimi model's param- 
eters, but does not contradict values of the parameter T^o ^ 3 • 10"'^'', i.e. excesses 
of an average energy of superthermal particles' component in 17(!) orders of 
equilibrium components' temperature. Corresponding to this parameter energy 
in the modern Universe at a temperature of the relic radiation 3°K is amount 
to 0, 8 * 10^GeV(!). Our next paper will be devoted to the investigation of the 
Universe model with the initial strongly violated thermal equilibrium. 

In conclusion authors express their thanks to professors V.N. Melnikov, Yu.S. 
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Vladimirov and D.V. Galtsov for a productive discussion of paper's matters. 
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